Disguising quantum channels by mixing and channel distance trade-off 
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We consider the reverse problem to the distinguishability of two quantum channels, which we call 
the disguising problem. Given two quantum channels, the goal here is to make the two channels 
identical by mixing with some other channels with minimal mixing probabilities. This quantifies 
how much one channel can disguise as the other. In addition, the possibility to trade off between the 
two mixing probabilities allows one channel to be more preserved (less mixed) at the expense of the 
other. We derive lower- and upper-bounds of the trade-off curve and apply them to a few example 
channels. Optimal trade-off is obtained in one example. We relate the disguising problem and 
the distinguishability problem by showing the the former can lower and upper bound the diamond 
norm. We also show that the disguising problem gives an upper bound on the key generation rate 
in quantum cryptography. 

PACS numbers: 03.67.-a, 03.67.Hk, 03.67.Dd, 03.67.Lx 
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I. INTRODUCTION 

Quantum information processing involves the transfor- 
mation of quantum states through quantum channels and 
it is often useful to quantify how far apart quantum states 
or quantum channels are. Depending on the problem 
at hand, different ways of measuring the distance may 
be adopted. Trace distance [l], and fidelity are 
two widely-used measures for quantum states. Trace dis- 
tance is particularly interesting because it corresponds 
to a measurement that distinguishes between two quan- 
tum states with the minimum error. Other distances for 
quantum states have also been studied recently, including 
the Monge distance @, the fcth operator norm Q, and 
the partitioned trace distance @|. For quantum chan- 
nels, measures Q have also been proposed based on ex- 
tending the fidelity measure [l(| and the trace distance 
measure fllj of quantum states. The diamond norm, in 
particular, is a trace-distance-based measure for quantum 
channels. It was first introduced in quantum information 
processing by Kitaev [ll[ for studying quantum error 
correction and has a nice operational meaning because 
it corresponds to minimum-error channel discrimination. 
As such, the diamond norm has been receiving a lot of 
attention since its introduction, in both the theoretical 
aspect fl2l - fl7j and the computational aspect [HI, [l9[ . 

While distinguishability (of quantum states and chan- 
nels) is a well studied problem, we consider the reverse 
problem - the disguising problem for quantum channels. 
Unlike the distinguishability problem in which the goal is 
to find a measurement that distinguishes between two (or 
more) states or channels, the aim in the disguising prob- 
lem is to find out the minimal mixing needed to make 
two (or more) quantum channels completely identical. In 
essence, this quantifies how much the effect of one chan- 
nel is partially carried out by another channel. In this 
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pilot study, we investigate the disguising problem for two 
channels. 

The diamond norm, which is a measure of distinguisha- 
bility, is not known to be easily computable, although a 
few algorithms based on semidefinite/convex program- 
ming and Monte Carlo strategies exist jT7il 22 . One mo- 
tivation for studying the disguising problem is to see if 
there is another measure for distances between quan- 
tum channels that is easily computable and conceptually 
similar to the distinguishability problem. Indeed, as we 
show in this paper, the disguising problem and the dis- 
tinguishability problem can be considered as dual to each 
other. We establish this by showing that the solution of 
the disguising problem can be used to lower and upper 
bound the diamond norm. This has an interesting impli- 
cation: the more distinguishable two quantum channels 
are, the more effort it takes to disguise one as the other. 
On the other hand, it turns out that solving the disguis- 
ing problem exactly is not easy either, but we are able to 
show how to efficiently compute lower and upper bounds 
of it. 

The disguising problem can be understood with the 
following operational interpretation. First note that the 
operational meaning of the diamond norm is based on 
the perspective of the receiver who tries to distinguish 
between two channels. A reverse perspective is to look 
at the channel intervener who tries to make the channels 
identical by minimal intervention. The channel inter- 
vener possesses the two original channels as black boxes. 
She is not allowed to open them and is only allowed to 
occasionally substitute each of them with some other ar- 
bitrary channel. We ask what are the minimal mixing 
probabilities needed to make the two intervened chan- 
nels identical? 

The precise problem statement is the following. Given 
two quantum channels £(p) = J^i EipE\ and J-(p) = 
J2i FipF} acting on an n-dimensional Hilbert space 
where Ei and Fi are n x n complex matrices represent- 
ing the Kraus operators of the channels, we consider the 
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processing 

£(p)^£'(p) = (l-p)£(p)+p£ A ( P ), (1) 
T(p) -> -F'(p) = (1 - + qT A (p) (2) 

such that £' = T' with the smallest p and q. In other 
words, we are interested in the least amount of substi- 
tution of £ and T needed to make them equal. This 
is illustrated in Fig. [TJ Operationally, the new channel 
probabilistically selects between the original channel and 
some other harmonizing channel, £a or Ta, which is yet 
to be determined. The smaller the mixing probability, p 
or q, the closer the new channel is to the original one. 
Thus, p and q serve as a distance between the two chan- 
nels. We note that in general the harmonizing channels 
£a and .Fa are not universal and depend on the original 
channels £ and T . In fact, the problem becomes trivial 
if we insist £a and T a to be universal for then £a = J~a 
(with cither p = q = 1 when £=/=J-orp = q = when 
£ = J-). Also, note that £' = J-' is trivially satisfied 
with £a = J~, J- a = £, and p + q = 1. This gives a linear 
trade-off between p and q. However, in general, better 
sub-linear trade-off can be obtained, as we show later. 
Note that £ a and J- a can be general quantum channels 
with arbitrary complexity. 

The pair of parameters (p, q) represents the trade-off 
between the two channels' mixing probabilities. The 
more mixing is imposed on one channel, the less mixing 
is required on the other. Note that if (p, q) is a trade-off 
point, (p + Sp, q + Sq) with Sp, Sq > is a also a trade-off 
point (and we call the former point strictly better than 
the latter). When given a region of achievable trade- 
off points, a trade-off curve can be obtained by tracing 
out the boundary such that no point is strictly better 
than another. This gives us a distance profile for the two 
channels. Thus, our measure is unique in that it is repre- 
sented by a 2-dimensional curve rather than a scalar as 
in other measures for quantum channels. On the other 
hand, a scalar distance may be obtained from our mea- 
sure in several ways, for example, (i) by imposing equal 
mixing probabilities p = q and regarding the minimum 
p = q as the distance between the two channels, or (ii) 
by regarding the minimum p + q as the distance. We will 
justify that these two are distances by showing that the 
triangle inequality holds. 

In this paper, we formulate the disguising problem as 
an optimization problem and show our main result in 
Sec. Mil Although solving it turns out to be difficult, we 
are able to obtain lower-bound and upper-bound on the 
(p,q) trade-off curve (cf. Eq. ([8])). We detail the anal- 
ysis of the bounds in Sec. HVl Next, we illustrate the 
computation of the bounds in a few examples for dif- 
ferent quantum channels in Sec. [V] In one special case, 
the analytical lower- and upper-bounds coincide, effec- 
tively producing the optimal trade-off curve. For the 
other cases, the numerically computed lower- and upper- 
bounds are quite tight, showing the effectiveness of the 
bounds. In Sec. IVIi we show that the disguising problem 
can lower and upper bound the diamond norm which 



quantifies the distinguishability of two quantum chan- 
nels. We also make three remarks about our distance 
profile for quantum channels in Sec. I VIII We discuss one 
application of the disguising problem in Sec. I Villi which 
is to bound the key generation rate in quantum cryptog- 
raphy. We conclude in Sec. IIXI 



II. NOTATIONS 

We denote a matrix A to be positive-semidefinite 
(PSD) by A y 0, the transpose of A by A 1 , and the 
conjugate transpose of A by A> . A is PSD if and only 
if A is Hermitian and its eigenvalues are non-negative. 
1 1 A 1 1 denotes the spectral norm of A which is the largest 
singular value of A or the largest eigenvalue of A if A is 
PSD. ||A||i = TtVaTA denotes the trace norm of A. 

B(T-L n ) denotes the set of all bounded linear operators 
in an n-dimensional Hilbcrt space TA. n . I p denotes the 
identity operator in a p-dimensional Hilbert space. A 
linear map £ : B(H n ) — s- B(H n ) is positive if £ (A) is 
PSD for all PSD A in B(H n ), and £ is completely positive 
(CP) if £ ® I p is positive for all positive integers p. £ is 
trace-preserving (TP) if Tr(£(A)) = Tr(A) for all A in 

B(H n ). 

A linear map £ : B(H n ) — > B(H n ) can be represented 
by a Choi matrix of size n 2 x n 2 [23( : 



C £ 



£(|0)(0|) f(|0)(l| 
£(|1)<0|) 

£(|n-'l><0|) ... 

n-l 

£ \i)(j\®£(\i)(j\). 

«',j=0 



£(|0><n-l|) 



£(|n-l>(n-l|) 



(3) 
(4) 



We define a function, which we call the channel sum func- 
tion, of the Choi matrix Cg of a linear map £ as follows: 



T{C £ ) := Tr*(C £ ) 
n-l 

= J2 \i)(j\-Tr[£(\j)(i 

i,j=0 



(5) 
(6) 



where Tr2 is the partial trace over the second system 
and t represents transpose. We remark that T(Cg) = 
I if and only if £ is trace-preserving (see Lemma [5] in 
Appendix Ia]) . 



III. PROBLEM FORMULATION AND MAIN 
RESULT 

A. Optimization problem formulation 

To solve for the optimal distance profile (p, q) of 
Eqs. ([lj and with the condition that the new channels 



FIG. 1. Two quantum channels, £ and T, are made identical by mixing the original channel £ (J-) with a harmonizing channel 
£a (J~a) with probability p (q). The designer is free to choose the non-shaded parts £a and J- a, but not £ and T. 



are identical, i.e., £' = J 7 ', we formulate the problem as 
minimize q (7) 
subject to £' = J 7 ', 



C £a h 0, 
C rA h 0, 
T(C fA ) = J ! 
T(C^ A ) = J ! 

where the minimization is overp, g, £a, and J-"a, for some 
fixed parameter /?. Here, we denote the Choi matrices of 
£a and Ja by Cg A and Cjr A , respectively (see Eq. ([3])), 
and T is the channel sum function defined in Eq. (|5j). 
The last four constraints demand that £a and J-"a be 
quantum channels (TPCP maps) (cf. Theorem [3] and 
Lemma [5] in Appendix [X| . Note that the roles of p and 
q are symmetric and can be interchanged when dealing 
with extreme points. Fig. [2] illustrates how to obtain the 
distance profile (p, q) by solving this optimization prob- 
lem with a range of /?. 

Note that for p = 1, the solution of q — is trivially 
obtained to make £' = T' . since we can choose £a = J~ 
in Eqs. ([T]) and ©. By the same token, (p, q) = (0, 1) is 
feasible. 

The distance profile (such as that in Fig. [2]) should be 
convex. This is because given two points (p, q) and {p 1 , q') 
that satisfies £' = T' [see Eqs. (Q} and ((5J)], any linear 
combination of them [i.e., ((1 — t)p + tp' , (I — t)g + tq') 
for some < t < 1] also satisfies it. It follows that any 
point on the line q — 1 — p is a feasible solution. 

B. Main result 

The solution of problem ((TJ) is indirectly bounded by a 
new variable as follows: 

n" 1 Tr[T(A+)] < a < min(||T(A+)|| , I). (8) 

Here, A + is a PSD matrix obtained by decomposing Cg — 
PCjr into the positive and negative subspaces by eigen- 
decomposition: 

C £ -PC T = A + -A_, (9) 
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FIG. 2. Distance profile obtained by solving problem (0. 
Among all the points (p, q) along the line corresponding to 
a fixed value of ft, we choose the point with the minimum q 
subject to the constraints of the problem. For each value of 
/3, we solve the optimization problem and obtain an optimal 
point (p, q). Repeating this for a range of j3 produces the solid 
curve. 



where A± are PSD matrices with support on orthogo- 
nal vector spaces (i.e., A+A_ = 0) and Cs and Cjr are 
the Choi matrices for the quantum channels £ and J 7 , 
respectively. 

The trade-off point (p, q) is obtained from a and /3 as 
follows: 

p=l-—L-, and (10) 

a + p 



Note that dp /da > and dq/da > implying that a 
smaller a gives rise to a "smaller" pair (p, q) in the 2- 
dimcnsional space. This means that the lower (upper) 
bounds of a in Eq. ([5]) obtained by varying (3 correspond 
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to a lower (upper) bound curve in the (p, q) space. 

Note that if the lower bound and upper bound of a 
coincide, the optimal a and thus optimal (p, q) are ob- 
tained. This happens if and only if A + already corre- 
sponds to a scaled quantum channel, i.e., T(A + ) = al, 
which is not the case in general. Also, note that Eq. (J5J) 
implies that the lower bound of a is always less than or 
equal to 1, and thus if T(A + ) = al, a < 1. 



IV. PROOF OF LOWER- 
UPPER-BOUNDS 



AND 



As noted earlier, the case of p = is trivial and thus 
we focus on the case p < 1 in the following. We now 
analyze problem (|7|), and as we will show later, directly 
solving this problem turns out to be difficult. Let us 
first focus on the condition that we want: £' = J-' . By 
Theorem [?] in Appendix fAJ we convert this condition to 
the Choi-matrix equivalence Cg> = Cjn, which implies 
that 



(1 - p)C e + pC £A = (l-q)C r + qCjr A 
1 „ V 



C £ - PCj 



l-p 



-C 



J~ A 



1-P 



-a 



(12) 
(13) 



where (3 = jztj;- We decompose the left-hand side into the 
positive and negative subspaces by eigen-decomposition: 

C £ -(3C^ = A+-A_, (14) 

where A± are positive semidefmite matrices with support 
on orthogonal vector spaces (i.e., A + A_ = 0). As such, 
by Theorem [3] in Appendix [Aj A± correspond to some 
CP maps. 

Note that G £ , Gj=, Gjr A , G £a , A + , and A_ are all 
Choi matrices. 

Comparing Eqs. (fT3|) and (fT4|) . since the positive and 
negative parts on the right-hand sides must match, the 
Choi matrices of the harmonizing channels must be of 
the form 



l-p 

P 
l-p 



Of a = A^ 
C £a = A_ 



X, 



X 



(15) 
(16) 



where X is some Hermitian matrix corresponding to the 
Choi matrix of some linear map. Note that A± may not 
correspond to scaled quantum channels because T(A±) 7^ 
al for some a > 0. The purpose of adding X is to make 
them scaled quantum channels so that T(A±+A) = a±I. 

Lemma 1. T(A+) = T(A_) + (1 - (3)1. 

Proof. Rearranging Eq. (|14|) and applying Corollary [21 
we have 

T(C £ + A_) = T(/3C^ + A+) 
=► T(G £ ) + T(A_) = T(J3Cf) + T(A+) 
=>-(l- J 8)/ + T(A_) = T(A + ) 



where we have used the fact that T(Cjr A ) = T(C £a ) = I 
since Cjr A and C £a are Choi matrices of quantum chan- 
nels. □ 

As a consequence, T(A + + X) = al if and only if 
T(A_ +X) = (a + /3— 1)1. Furthermore, from Eq. (fT5|) . 
since T(C^- A ) = /, we have 



(17) 



l-p 



The same expression is obtained when we consider 
Eq. (|16p with T(C £a ) = I. Thus, minimizing q given 
/3 fixed is equivalent to minimizing a given (3 fixed, since 



a 



1 



(18) 



is an increasing function of q. The original problem ([7|) 
becomes 

a — minimize a (19) 
subject to A + + X y 0, 

a_ + x y 0, 

T(A + + X) = al, 

T(A_ + X) = (a + - 1)1, 

where the minimization is over Hermitian matrix X given 
(3 fixed, and A± are from Eq. (fT4")) . Note that the fourth 
constraint is redundant due to Lemma [1] and is shown 
only for completeness. Once a is found, we can compute 
p and q from Eqs. (JT7|) and ([15]). 

We investigate the form of X. Since Cjf a and C £a 
represent quantum channels, they are PSD. This means 
that, according to Eqs. Jl5|) and fT6]). A± + X are PSD. 
However, this does not mean that X is also PSD, and this 
makes finding the optimal X difficult. Nevertheless, we 
have the following constraint on X which helps us bound 



Lemma 2. 

Tr(X) > 0. 



The constraints of problem (|19p implies 



Proof. Since A + and A_ are the positive and negative 
subspaces of the matrix of Eq. (fT4]) . we can identify 
non-overlapping projectors P+ and P_ onto them respec- 
tively. We also define the projector onto the remaining 



subspace Pq 
have 1 



■■I-P+-P-. Since A± + X is PSD, we 

Tr[P_(A + +A)] >0, 
Tr[P+(A_ + X)] > 0, and 
Tr[P (A±+A)] >0, 



which implies that 

Tr(P_A) > 0, 
Tr(P+A) > 0, and 
Tr(PoA) > 0. 

Summing these terms gives the desired result. □ 
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This lemma implies that the non-zero eigenvalues of 
X cannot be all negative, but X can have positive and 
negative eigenvalues. 



Theorem 1. The optimal value of problem ([19} is upper 
bounded by ||T(A+)||. 

Proof. To show an upper bound, we only need to find 
a feasible X. Choose M = ||T(A+)p - T(A+) as the 
difference between two channel sums. Certainly, M is 
PSD and thus can be written as M = D\d$ where Dq 
is a square matrix. M represents the channel sum of 
the channel p — > DqpDq. Let X = \Do)(Dq\ be the 
Choi representation of this channel where \Do) is the 
vector form of D (cf. Eq. {53}). Since X is PSD, 
A± + X is PSD and the first two constraints of prob- 
lem (Tr5} are satisfied. Note that T(X) = M by con- 
struction (cf. Lemma U and Definition [1}. Therefore, 
T(A+ + X) = T(A+) + T{X) = ||T(A+)||J by Corol- 
lary [21 

Note that for this upper bound, we have chosen X to 
be PSD. □ 

We computed ||T(A + )|| for random quantum channels 
and found cases with ||T(A + )|| > 1. Nevertheless, a < 1 
is also a valid bound. 



Lemma 3. The optimal value of problem (|19|) is upper 
bounded by unity, i.e., a < 1. 

Proof. Set the harmonizing channels in Eqs. ([T}-([2} to be 
£a(p) = HP) and F A {p) =£{p). Then, £'(p) = P(p) 
is satisfied with q = 1 — p, which means that a = 1 
according to Eq. (|TT} . Based on Eq. (fl4|) . we set X = 
C £ -A + = /30f-A_. Then, we have A+ + X = C £ h 
and A_ + X = PCjr y 0. As such, the constraints of 
problem (1191) are satisfied with a = 1. □ 



We remark that in the above proofs of the two upper 
bounds, we have explicitly constructed X. Therefore, 
problem (|19p is always feasible. 



Theorem 2. The optimal value of problem (TT5} is lower 
bounded by n _1 Tr[T(A + )]. 

Proof. The channel sum is T(A + + X), and the sum its 
eigenvalues is 

Tr[T(A+ + X)} = Tr[T(A + )] + Tr[T(X)] 
> Tr[T(A + )], 

where the first line is due to linearity of T (cf. Corol- 
lary [3} and the second line is due to Corollary [1] and 
Lemma H which imply Tr[T(JQ] = Tr(X) > 0. Finally, 
since T(A + + X) = al (cf. problem |19}), we have 
a > n- 1 Tr[T(A + )]. 

□ 

In summary, the solution of problem fj 19[) is bounded 
as follows: 



If A_)_ already corresponds to a scaled quantum chan- 
nel, i.e., T(A+) = al for some a, then the optimal solu- 
tion can be found: a = 77,- 1 Tr[T(A + )] = ||T(A+)|| = a. 
In this case, Cjr A = a -1 A + and Cg A = (a + $ — 1) _1 A_ 
can be found from Eqs. (fTS} and (fT6|) with X = and 
Eq. (TTT}. 



A. Procedure for computing the lower- and 
upper-bound (p, q) curves 

Suppose that we are given two quantum channels £ 
and T of dimension n. 

1. Compute the Choi matrices Cs and Cjr for the two 
channels using Eq. ([3}. 

2. Fix P in the range of (0,oo). [Note that p = or 
P = oo corresponds to q = 1 or p = 1 respectively, 
and these are trivial cases because either £' or T' 
becomes arbitrary] 

3. Eigen-decompose Eq. (TTJJ to obtain A±. 

4. Compute the channel sum T(A+) using Eq. @. 

5. Compute the lower and upper bounds on a using 
Eq. (HQ}. 

6. Given a bound, denoted as a, solve for p and q 
using Eqs. JTO} and (|TT) . 

We can repeat this procedure for a range of p to obtain 
the lower- and upper-bound (p, q) trade-off curves. 



V. EXAMPLES 

A. Difference between bit-flip and phase-flip 
channels 

Given the bit-flip and phase-flip channels, 

£{ P ) = {l-a)hph + aX P X (21) 
F{p) = (1 - b)I 2 ph + bZpZ (22) 

where a and b are the bit-flip and phase-flip probabilities, 
and 



"1 0" 


, x = 


'o r 




1 " 


, z = 


1 


1 




-1 



we compute the Choi matrices for the two channels and 
find the difference 



C £ - PC T 

=(l-a-p + bp)\e 1 )(e 1 \ - bp\e 2 )(e 2 \ + a\e 3 )(e 3 \ 



(23) 



n" 1 Tr[T(A + )] < a < min(||T(A+)|| , 1) 



(20) 



where (e x | = [1,0,0,1], (e 2 | = [1,0,0,-1], and (e 3 | = 
[0, 1, 1,0]. Next, we separate this into the positive and 
negative subspaces as in Eq. {H}. Note that since we 
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consider a,b,/3 > 0, the second term of the equation is 
negative and the third term is positive, while the first 
term can be non-negative or negative. 

Case 1: 1 - a - (3 + b(3 > 0. According to Eq. (fig]), we 
have 

A+ = (1 - a -p + bp)\e x )(ex\ + a|e 3 )(e 3 |, 
A_ = &/?|e 2 )<e 2 | 

and 

T(A+) = (1-/3 + 6/3) J 2 , 
T(A_) = &£I 2 . 

Therefore, using the bounds in Eq. ([20]) . we obtain the 
optimal solution of problem (fT9|) as a = 1 — f3 + b(3. 

Case 2: l-a-/? + &/?<0. According to Eq. ([13]). we 
have 

A + = a|e 3 )(e 3 |, 

A_ = -(l-a-/3 + 6 ( S)|ei)(e 1 | + 6/3|e 2 )(e 2 | 

and 

T(A+) =aJ 2 , 

T(A_) = (-l + a + /3)7 2 . 

Therefore, using the bounds in Eq. ([20]) . we obtain the 
optimal solution of problem (|19[) as d = a. Note that 
we are able to obtain the optimal solution in both cases 
instead of upper and lower bounds. 

Finally, with a found for each case, we can compute a 
relation for p and q using Eqs. (|17[) and (TT5]> : 

(p = b-bq , ifl-a-0 + bj3>O . . 

\q = a-ap , ifl-a-j9 + 6j8<0' 1 j 

This relation is depicted as the solid curve in Fig. [3] where 
the top-left (bottom-right) part corresponds to the first 
(second) case in Eq. ([24]) . Essentially, the cusp in the 
figure is due to the transition from case 1 with 2 positive 
and 1 negative eigenvalues to case 2 with 1 positive and 
2 negative eigenvalues in Eq. ([23]) . 

B. A pair of random qubit channels 

We randomly generated two qubits channels each hav- 
ing four Kraus operators and they are listed in Ap- 
pendix [B] Using the procedure given in Sec. IIV Al we 
compute the lower- and upper-bound curves which are 
shown in Fig. [4] We make two observations. First, four 
cusps are obvious in the lower-bound curve, which are 
due to the transition of an eigenvalue of Eq. ([14]) from 
positive to negative (or vice versa). Note that at most 
four cusps can occur since the dimension of the Choi ma- 
trices are four. Second, there are regions where q = for 
a range of p and where p = for a range of q (the former is 
much bigger than the latter). These regions correspond 
to the case that one channel contains another channel 
and we will clarify this concept in Sec. IVIIBl later. 
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FIG. 3. Optimal trade-off curve for the mixing probabilities 
p and q [defined in Eqs. |T]) and @] for the bit- and phase- 
flip channels given in Eqs. (|21|) and (|22[) . This curve is the 
solution to problem ([7]) or problem (|19p . 




FIG. 4. Lower- (purple) and upper-bound (blue) curves for 
the mixing probabilities p and q [defined in Eqs. {TJ and ((2|] 
for two random qubit channels. 



C. A pair of random four-dimensional channels 

We randomly generated two four-dimensional channels 
each having four Kraus operators. (We do not list them 
here as they take up a lot of space.) Fig. [5] shows the 
bounds. We remark that the solid upper-bound curve 
computed with Eq. (|20|) is not useful since it is above the 
q = 1 — p line which is a trivial upper-bound. Neverthe- 
less, the lower-bound curve is useful since it allows only 
a narrow gap with the upper-bound line. 

This example brings up an important point: in general, 
we should take the convex hull of an upper-bound curve 
as a refined upper-bound. 

We also remark that the maximum number of cusps 
in the bounding curves is 16 since the dimension of the 
channels is 4. However, they are not apparent in the 
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FIG. 5. Lower- (purple) and upper-bound (blue) curves for 
the mixing probabilities p and q [defined in Eqs. |T} and @] 
for two random four-dimensional channels. The solid curves 
are the bounds computed with Eq. (|20p and the dashed curve 
is q = 1 — p. 



figure. 

VI. RELATION WITH DISTINGUISH ABILITY 

The diamond norm is related to the minimum error 
in discriminating between two quantum channels and is 
defined as 

||S || ^ maX || (/£«>£)(/>) ||i. 
p 

Here, an ancillary Hilbert space K, is introduced and Ik. 
is the identity map acting on it. The dimension of lie, is 
the same as the dimension of the Hilbert space of Q [Tl| . 
The minimum error in distinguishing between £ and T 
is given by [24[ 

P(error) = i(l-||f-^y. 



A. Upper bound 

We can upper bound the diamond norm with the mix- 
ing probabilities p and q of our disguising problem as 
follows: 

\\£-r\\<> = ll(i - p)£ - (i - 1)? + p£ - iHo 

— \\q^A — p£a + p£ - qF\\o 

< gll^Allo +p\\£a\\o +p||£||<> + gll^llo 

= 2(p + q), (25) 

where the second line is due to Eqs. (QJ and ([2]) with the 
disguising condition £' = T' satisfied, the third line is 
due to the triangle inequality of the trace distance, and 
the fourth line comes from the fact that the trace norm of 



the channel output (a density matrix) is one. Note that 
Eq. (|25|) holds for any feasible (p,q) satisfying £' = J 7 ', 
not just the optimal (p, q) trade-off curve. 

When the two channels £ and T are perfectly distin- 
guishable, || £ — J 7 ^ = 2. On the other hand, in our 
disguising problem, p + q = 1 is always achievable in 
Eqs. (p} and @ since we can set p = 1 — q, £a = J 7 , and 
.7-a = £■ Therefore, Eq. is tight in this case. 



B. Lower bound 

We focus on the case where p = q, which means that 
(3 = 1. In this case, from Eq. flSJ), the smallest p = q 
satisfies 



P 



1 



P 



1 



<I|T(A + )||, 



(26) 



where A + is the positive subspace of Cg — Cjr. 

We divide A + (of size n 2 xn 2 ) into n x n blocks of equal 
size and denote the block as A+y. Thus, A + = 

Y^ijlo K)0'l ® A+ij, and it follows from the definition of 
T in Eq. © that 

l|T(A + )|| = 



n-l 



£ |i)(j|-Tr[A +ii ] 
»>j=o 

n— 1 n — 1 

max (4>\i)(i\<t>) ■ (zk\&+ij\zk) 
\4>) i,j=o fe=o 
n-l 

max (0| ® (z k \A + \<f>) ® \z k ) 



k=0 



< n max 

W.I*} 



l*>> 



(27) 



where on the second and third lines {|z/c)} is an orthonor- 
mal basis, and (<f>\<f>) = (z\z) = 1. 

Next, we consider the diamond norm: 



>max||(/®(£:-^))(|^)^|)|| 

m 

= max \{a\(I®(£ -7)) (ty) (VI) k) 

|cr)#) 



(28) 



where the second line is due to the fact that the trace 
norm is no less than the spectral norm, and both the 
auxiliary system and the original system have dimension 
n. Without loss of generality, using the Schmidt decom- 
position, we can express 



|<7) = ^7«l^)h) 
i=0 
n-l 



i=0 
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where {|-Bi)} is an ortho normal basis, 



equal to 



(cj|o"i) = (ipi\ipi) — 1 for all i = 0, . . . , n — 1, and 

n — 1 n — 1 

Ew 2 = En 2 = 1 - 



(29) 



n-l 



i=0 i=0 

Continuing with Eq. 



the term to be maximized is where 
I 



K«|G|«)| 



G 



(£-T)(\iPi)(iPo\) 



(£- J-)(|V„-i)(Vn-i|). 



and 



ko>A£7o 



|cr„_i)A*_i7n-i 



Note that G is not a standard Choi matrix since there is 
no requirement that {|?/>i)} is an orthonormal basis and 
also (v\v) = J2i |A;| 2 |7i| 2 may not be unity. 
Continuing with Eq. (|28p . we have 



Q = max max |(i>|G|v)| 

{ki)}>{7i} {\1k)},{>«} 

subject to the constraints in Eq. (|29p . Since \ipi) = \i) 
and A; = l/V"- satisfy the constraints, 



C. Summary 

Using Eqs. (|25p and (|3ip . the smallest p 
satisfy 



V 



n 2 (l-p) 



< \\S-T\L<4p. 



q must 



(32) 



This shows that the disguising problem and the distin- 
guishability problem are dual: when two channels are 
easy to distinguish (the diamond norm is large), it re- 
quires great effort to disguise one channel as the other 
(p = q is large); and the reverse also holds. Note that 
the lower bound of Eq. (|32p is less than or equal to the up- 
per bound since p = q is at most 1/2 due to the fact that 
any point on the line q = 1 — p is feasible (cf. Sec. lIII Al . 



Q>max-\{v\[C e -Cjr)\v)\ 
\v) n 

= l -\\c £ -c A 

= imax(||A + ||,||A_||), 

where the maximization is over any vector \v) with 
(v\v) = 1 and G with the substitution = \i) is equal 
to Gf — Cjr. Finally, note that ||A+|| is larger than or 
equal to the maximization term in Eq. (|27p . Therefore, 

l|£-^IU> Al|T(A + )|| (30) 

and combining with Eq. (|26p . the smallest p = q must 
satisfy 



VII. OTHER REMARKS 

A. Triangle inequality 

We apply the notion of triangle inequality to our mix- 
ing probabilities. Suppose £ and T are compatible with 
mixing probabilities (p, q) and Q and T with (p',q'), 
meaning that 

(1 - p)£(p) + p£ A (p) = (1 - q)T(p) + qJ- a(p), (33) 

(i -p')g(p) +p'Ga(p) = (i - q')Hp) + q'^kip)- (34) 

Note that the harmonizing channels J-"a and T' A are dif- 
ferent in general. We want to infer the distance profiles 
{p",q") for £ and Q from the distance profiles (p 7 q) and 
(p' 7 q'). To do this, we propose the following method: 
cross-multiply Eqs. (|33[) and (|34[) to make the coefficients 
of J- equal and add additional terms to the two resultant 
equations to make the overall harmonizing channels on 
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the right-hand sides equal. The result is 

(1 - q') [(1 - p)£ + p£ A ] + (1 - q)q'F' A 

= (1 - q') [(1 - q)T + qT A \ + (1 - q)q'F' A , and 
(1 - q) [(1 - pOS + j/0 A ] + (1 - g')g^A 

= (1 - g) [(l - ^ + q F' A \ + (l - 9 ')^a, 

where we drop the dependence on p for simpler notation 
and assume not both q and q' equal to 1. Thus, the two 
left-hand sides are equal, giving 



(1 - q') [(1 - p)£ + p£ A ] + (1 - q)q'J" A 

= (i-q)[(i-p')g +p 'g A ] + (i 



q')qT A . 



This equation is interpreted as £ occurring with proba- 
bility (1 — q')(l—p) /(1 — qq') and its harmonizing channel 
with probability 



„ = p(i - q') + (l - q)q' 



l-qq' 

and C/ occurring with probability (1 — g)(l —p')/(l 
and its harmonizing channel with probability 

„_ p'(l- g ) + (l-g')q 
l-qq' 



(35) 



qq 



(36) 



With Eqs. (|3"5j) - (|3"rJ)) . given an achievable pair of mixing 
probabilities (p, q) for £ and T and another pair (p',q') 
for J 7 and G, we can compute an achievable pair (p", g") 
for £ and Tracing out the entire distance profiles (p, q) 
and (p',q') produces an achievable region (p",q"). Note 
that this region is not in general a curve. Nevertheless, 
the bounding curve to this region can be regarded as 
a distance profile for £ and Q. This curve is certainly 
achievable but may not be optimal and so it represents 
an upper bound to the optimal trade-off distance profile. 

Note that even though we have the assumption qq' < 1 , 
we can still obtain the end points. In particular, when 
(p,q) = (0,1) and q' < 1, we get (p",q") = (0,1); when 
( P ',q') = (0,1) and q < 1, we get (p",q") = (1,0). 

As an example, consider £ and T given in Eqs. (|21|l 
and (HH and 

Q{p) - (1 - c)I 2 ph + c(XZ)p{ZX) 

with a = b = c = 0.2. Figure [6] shows the achievable re- 
gion for £ and Q obtained by Eqs. (|35|) and (j36j) together 
with the optimal curve obtained in Sec. IV Al Note that 
due to symmetry, any pair of £, J-, and Q has the same 
optimal trade-off curve. 

We now consider the triangle inequality for two scalar 
distances derived from our 2-dimensional measure: (i) 
the minimum of p subject to p = q and (ii) the minimum 
of p + q. For case (i), suppose that the minimum for £ 
and T is po with go = Po and the minimum for T and 
G is p' with q' = p'q. Substituting these four parameters 
into Eqs. (|35|) and (f36j) gives 



= q 



P = 



Po ( 1 ~Po) + (1 - Po )p'o 
1 - PoPo 




FIG. 6. Achievable region for £ and Q obtained by Eqs. p5[) 
and (|36p using the mixing probabilities of £ and T and the 
mixing probabilities of J- and Q. The dashed curve (red) is 
the optimal trade-off. 



This means that the p = q condition is preserved. Since 

(1 -p'o)i l - PoPo) < 1 and i 1 -Po)(l - PoPo) < 1. we 
have p" < po +Po- Finally, since p" is only an achievable 
upper bound to the optimal mixing probability for £ and 
Q, the triangle inequality is satisfied. 

For case (ii), suppose that the minimum for £ and T 
is po + go an d the minimum for T and Q is p' + q' . 



Substituting these four parameters into Eqs. Q35[) 
([36]) gives 

(Po + 9o)(l-g ) + (Po + <Zo)(l-9o) 



and 



P 



q" = 



qoq 



Using the same argument as in case (i), we have p" + q" < 
(Po + '?o) + (Po"'"9o) an d t ne triangle inequality is satisfied. 



B. Channel containment 

Given two quantum channels £ and J 7 , we introduce 
the notion that £ contains J 7 if J 7 is part of a mixture of 
£: 



£{p) = (1 - q)F(p) + qF A (p) 



(37) 



where J- A is some quantum channel and wc require that 
q < 1 to make the containment of T non-trivial. To see if 
Eq. (|37p holds, we convert it to the Choi representation 
and proceed as in Eq. (fT3|) with p = 0. The question 
becomes for what values of q is C$ — (1 — q)Cjr = gCjF A 
PSD. This is the same as asking whether A_ = in 
Eq. (|14[) . In general, we fix a value of q and perform 
eigen-decomposition to see if A_ = 0. But for the 
case that Cjr is invertible, we can find the minimum 
q by equating 1 — q with the minimum eigenvalue of 

c~ 1/2 c £ c~ 1/2 . 

Note that Cjt a is automatically trace-preserving, since 
gT(C^ A ) - T(C £ - (1 - g)C7^ = I - (1 - q)I = ql as 
both £ and J- are trace-preserving. 
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C. Composition of quantum channels 

Suppose that Ei and Ti can be made compatible with 
Pi and qi according to the processing in Eqs. (P)-©, 
where i = 1,2. This means that £[{p) = T-(p) or 

(1 - Pi)Si(p) +Pi£Ai{p) = (1 - qi)Fi{p) + Qi^Aiip) 

for all density matrices p and i = 1,2. Then the 
composed quantum channels £ 2 ° £\ and T 2 o T\ can 
also be made compatible with p = p\ + P2 — P1P2 and 
q = qi + q 2 - qiq 2 : 

(1 - p)£ 2 o fi(p) + P £a(p) = (1 - q)T 2 a Ti{p) + qT A (p). 

(38) 

To show this, note that since £i and Ti are density 
matrices, the following holds: 

£ 2 (£[(p))=T 2 (T[(p)). (39) 

Expansion of the LHS gives 

£' 2 o£[ ={l-p 1 ){l-p 2 )£ 2 o£ x 

+ Pi(l - P2 )£2 o £ A i (40) 

+ (1 -pi)p 2 ^A2 O £l 

+ PiP 2 £a 2 ° £ai- 

We can readily see that the first term is the origi- 
nal composed channel with mixing probability 1 — p = 
(1— pi)(l— p 2 ) and the sum of the last three terms repre- 
sents the harmonizing channel £a with mixing probabil- 
ity p — pi+p 2 — Pip 2 - Together with a similar argument 
for T'2 o T[ proves the claim. 

We can also argue that the composed quantum chan- 
nels £2 o £ 1 and T 2 T\ can be made compatible with 
p = Pi+p 2 and q = q\ + q 2 , by breaking up the first term 
of Eq. (|40|) and allocating the portion p\p 2 £ 2 o £\ to the 
harmonizing channel (and similarly for T). 



VIII. APPLICATION TO QUANTUM 
CRYPTOGRAPHY 

The disguising condition £' = T' can be used to upper 
bound the key generation rate in quantum cryptogra- 
phy [25|, [2(| . The intuitive idea is that the more easily 
the eavesdropper's channel can be disguised as the legit- 
imate user's channel, the smaller is the amount of the 
generated key. We establish this idea quantitatively re- 
lating the mixing probability and the key generation rate. 

In quantum cryptography with one-way forward com- 
munications, Alice repeatedly sends a quantum state p a 
to Franky to establish a secret key, where a = (ao,ai) 
and ao (a\) is Alice's raw key basis (value) chosen in- 
dependently between different transmissions. After the 
reception of the sequence of states by Franky, Alice sends 
classical information (including basis information, error 



correction information, and privacy amplification infor- 
mation) to him in order to correct bit errors and remove 
any information the eavesdropper Eve may have on the 
final key . Suppose that Eve launches a collective at- 
tack [27| which means that she applies the same unitary 
transformation U to each state sent by Alice (with suffi- 
cient ancillas). Thus, Franky's channel and Eve's channel 
are given as follows: 

T(p a ) = Tr E (U(p a ®\0}{0\)U^ (41) 
£{p a ) = Tr F (J7(/o a ® \Q){Q\)U f ) (42) 

where we assume without loss of generality (w.l.o.g.) that 
the entire Hilbert space is divided into two systems E and 
F. Furthermore, for simplicity and w.l.o.g., we assume 
that E and F have the same dimensions n (which can be 
assured by padding zeros as needed). 

A key rate upper bound is given by the classical secret 
key capacity formula [28| 



R = sup I(U;F\V)- I{U-E\V). 



(43) 



Note that the use of the classical formula is valid since 
when considering the upper bound, we can assume one 
particular strategy of Eve, which is to measure her quan- 
tum states separately. Here, A is Alice's random variable 
holding the raw key a, F and E are the Franky's and 
Eve's random variables holding the measurement out- 
comes / and e respectively. 

We consider the upper bound of the key rate for the 
case where the disguising condition is 



T(p) = (1-p)£(p)+ P £a(p). 



(44) 



That is, 1—p fraction of Franky's channel is Eve's chan- 
nel. Thus, one would expect that only the remaining 
fraction of p could be used to generate a secret key and 
the key rate would be on the order of p. 

For simplicity of discussion, instead of Eq. (|43j) . we 
bound the key rate expression without the processing: 



R' = I(A; F)-I(A;E). 



(45) 



Suppose that Bob's POVM is {Mf,a } where 
Mf_ ao = I and it is dependent on the raw key basis 
a . Applying this POVM to Eq. (|44l) produces classical 
probability distributions 

P AF {aJ) = PA(a)Tr(M /iQ0 J-(p a )) 
P AE {aJ) = P A {a)Tr(M f , ao £(p a )) 
PaeAvJ) = P A (a)Tr(M f , ao £A(pa)) 
which are related by 

P AF (aJ) = (l-p)P AE (aJ)+ P P AEA (aJ). (46) 

Note that this relation can be explained by the follow- 
ing hypothetical probability distribution 



P A Fz(a,f,z) 



(l-p)P AB (a,f), ifz = 
pP AE& (a,f)i if 2 = 1 
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in that J2z=o 1 Pafz(o>, f, z) is equal to Eq. (|46l) . 
Now, we bound the first term of Eq. ([45]) as follows: 



I(A;F)<I(A;FZ) 

= I(A;F\Z)+I(A;Z) 

= (l-p)I(A;F\z = 0)+ P I(A;F\z = l) 

= (l-p)I(A;E)+pI(A;E A ) 

where on the second line we have I{A; Z) = since 
P AZ {a,z) = P A (a)P z (z). Therefore, using Eq. (g5]), the 
key rate is bounded as 

R' < (1 - P)I{A; E) + pI(A; E A ) - I(A; E) 
= p[I(A;E A )-I(A;E)} 

<p\og 2 n (47) 

where the last inequality is due to Hole vo- Schumacher- 
Westmoreland channel capacity theorem [29|, |3(| and the 
fact that the maximum entropy for an n-dimensional 
state is log 2 n. A similar analysis can be applied to the 
original key rate expression in Eq. t)43[) to obtain the same 
upper bound in Eq. (|47|) . 



IX. CONCLUSIONS 

The disguising problem tries to make two quantum 
channels identical, which is the reverse of the distin- 
guishability problem which tries to maximize their dif- 
ference in the measurement statistics in order to discrim- 
inate between them. Indeed, we showed that the two 
problems are related by proving that a certain combina- 
tion of the mixing probabilities of the disguising problem 
upper bounds the diamond norm of the distinguishabil- 
ity problem. We also showed that the triangle equality 
holds for two scalar distances derived from the mixing 
probabilities. 

Conventional measures on quantum channels are 
mostly based on trace distance or fidelity which are both 
concepts derived from measuring the distance between 
quantum states. In this paper, we propose a new measure 
genuinely for quantum channels. Note that the applica- 
tion of our measure to quantum states is possible but the 
result would be rather trivial since there is no more the 
need of making sure the harmonizing channels satisfy the 
TP condition for a linear map (which is ensured by the 
addition of X in Eqs. (fT5|) and ([To]) ). This extra condi- 
tion makes the calculation of our measure for quantum 
channels difficult. Nevertheless, we obtain lower- and 
upper-bounds curves which are close in many cases. 

Our measure is based on the notion of minimizing the 
probabilities of channel mixing, which can be viewed as 
the costs for an channel intervener to make two channels 
the same. We show how these costs are linked to the key 
generation rate in quantum key distribution. The inves- 
tigation of how these costs are linked to other quantum 
information processing tasks is a topic for future research. 



Also, open problems include efficient /approximate com- 
putation of our measure, and the effect of extending the 
Hilbert space dimensions of the channels by including an- 
cillary systems (i.e., £ becomes £ <E> I) in the calculation 
of our measure. 
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Appendix A: Useful results related to quantum 
channels 

This section discusses the tools and definitions re- 
lated to quantum channels. A quantum channel is a 
linear map that is completely-positive (CP) and trace- 
preserving (TP). 

Theorem 3. (Choi's theorem (23[) Given a linear map 
£ and its Choi matrix Cg, Cg is PSD if and only if £ is 
a completely-positive map. 

When Cg is Hermitian, £ can also be represented in 
the operator-sum form: 



£(p) = Y,\iE iP E}, 



(Al) 



where Xi £ K and Ei £ C n ' n . This can be seen by taking 
some decomposition (e.g., eigen-decomposition) of Cg to 
be 



C £ =Y J h\E i )(E i 



(A2) 



where A; 6 1 and |Bj) € C" and rearranging the 
vector \Ei) into the square matrix Ei as follows: 

E t (2,l) 

E t {n,l) 

Ei(l,2) =vec ( Si ) (A3) 
Ei{n,2) 
Ei(n,n). 

where Ei(k,l) is the (k,l) entry of Ei, and the vec op- 
erator creates a vector by stacking the columns of its 
operand (see, e.g., Ref. |3l|)- The dimension of Ei is 
n x n and the dimension of Cf is n 2 x n 2 . Note that the 
operator-sum form is not unique; there can be more than 
one such form corresponding to the same Choi matrix. 
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Observation 1. The Choi matrix of any channel given 
in the operator-sum form of Eq. (|A1|) can be constructed 
by using Eq. Q or Eq. (|A2[) . The equivalence of these 
two ways can be checked easily by direct expansion. 

The next Theorem follows directly from the definition 
of the Choi matrix in Eq. ([3]) . 

Theorem 4. £{p) = F(p) for all density matrices p if 
and only if C £ = Cj=. 

Definition 1. We define the channel sum for linear map 
£ with a Hermitian Choi matrix as 



(A4) 



which has dimension n x n. 



The next lemma shows how to obtain the channel sum 
from channel outputs directly. 

Lemma 4. 

n-l 

Ts=J2 K>(i|-Tr[£(|j)(i|)] 

i,j=0 

= Tr*(C £ ) 

where Tr2 is the partial trace over the second system and 
t represents transpose. Note that here, it does not matter 
whether the transpose is taken after or before the partial 
trace. 

Proof. The element of S £ defined in Eq. (|A4[) is 
(i\T £ \j) = {i\J2*kEtE k \j) 



Tr 



J2*kE k \j){i\El 



T±[£(\j)(i\)] 



□ 



Therefore, T £ is independent of the operator-sum form 
of £ and is dependent only on the Choi matrix. This 
allows us to define the channel sum function 



T(Cg) :=T4(C £ ) 



(A5) 



where we used Lemma |4] The introduction of T facili- 
tates the discussion of the channel sum with reference to 
only the Choi matrix. 



Lemma 5. A linear map £ is trace-preserving if and 
only if T(C £ ) = I. 

Proof. For £ to be trace-preserving, the following must 
hold for all density matrices p: 

Tr(p) = Tr[£(p)} 
= Tr(T£ P ) 

where the last equality is due to Eq. (|A4[) . Since this 
holds for all p, T £ = I. Then, using Eq. (|A"5|) . T(C £ ) = I. 
The proof for other direction is obvious. □ 

We remark that the eigenvalues of a Choi matrix C 
and of its channel sum T(C) are in general not the same. 
Nevertheless, they have the same trace. 

Corollary 1. Tr(C) = Tr(T(C)) for some Choi matrix 
C. 

Thus, the trace of the Choi matrix of a quantum chan- 
nel is its Hilbert space dimension n because the channel 
sum of a quantum channel is /„. This corollary will be 
useful when we consider the lower bound of the channel 
distance. 

Corollary 2. C £ +C? = C £+T and T{C £+T ) = T{C £ ) + 
T(Cjr) for linear maps £ and J 7 . 

Also, note that it can easily be checked that if C is 
PSD, T(C) is also PSD. 



Appendix B: Channel specification for Sec. IV Bl 

The two channels are £{p) = Xa=i EipEj and F{p) 
Ei=i F iP F i: where 



E 1 = 

E 2 = 

E 3 = 

E± = 

Fx = 

F 2 = 

F 3 - 

Fi = 



-0.504828 
-0.0133105 



-0.331944 
0.295026 



0.419485 0.158018 
0.330761 0.0616354 

0.464696 0.251826 
-0.312786 0.165248 

0.160149 -0.346665 
-0.346665 0.750403 

-0.20917 -0.248828 
0.382771 -0.451866 

-0.62412 -0.425856 
0.286902 -0.0613943 

0.216184 -0.422341 
-0.403389 0.451605 

0.236514 0.269256 
0.269256 0.306531 
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